Computational basis of SARS-CoV 2 main protease inhibition: an insight from molecular dynamics simulation based findings.

The coronavirus disease 2019 (COVID-19) pandemic is caused by newly discovered severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2). One of the striking targets amongst all the proteins in coronavirus is the main protease (Mpro), as it plays vital biological roles in replication and maturation of the virus, and hence the potential target. The aim of this study is to repurpose the Food and Drug Administration (FDA) approved molecules via computer-aided drug designing against Mpro (PDB ID: 6Y2F) of SARS CoV-2 due to its high x-ray resolution of 1.95 Å as compared to other published Mprostructures. High Through Virtual Screening (HTVS) of 2456 FDA approved drugs using structure-based docking were analyzed. Molecular Dynamics simulations were performed to check the overall structural stability (RMSD), Cα fluctuations (RMSF) and protein-ligand interactions. Further, trajectory analysis was performed to assess the binding quality by exploiting the protein-residue motion cross correlation (DCCM) and binding free energy (MM/GBSA). Tenofovir, an antiretroviral for HIV-proteases and Terlipressin, a vasoconstrictor show stable RMSD, RMSF, better MM/GBSA with good cross correlation as compared to the Apo and O6K. Moreover, the results show concurrence with Nelfinavir, Lopinavir and Ritonavir which have shown significant inhibition in in vitro studies. Therefore, we conclude that Tenofovir and Terlipresssin might also show protease inhibition but are still open to clinical validation in case of SARS-CoV 2 treatment.Communicated by Ramaswamy H. Sarma.

Glycoprotein molecular dynamics analysis: SARS-CoV-2 spike glycoprotein case study.

Molecular Dynamics (MD) is a method used to calculate the movement of atoms and molecules broadly applied to several aspects of science. It involves computational simulation, which makes it, at first glance, not easily accessible. The rise of several automated tools to perform molecular simulations has allowed researchers to navigate through the various steps of MD. This enables to elucidate structural properties of proteins that could not be analyzed otherwise, such as the impact of glycosylation. Glycosylation dictates the physicochemical and biological properties of a protein modulating its solubility, stability, resistance to proteolysis, interaction partners, enzymatic activity, binding and recognition. Given the high conformational and compositional diversity of the glycan chains, assessing their influence on the protein structure is challenging using conventional analytical techniques. In this manuscript, we present a step-by-step workflow to build and perform MD analysis of glycoproteins focusing on the SPIKE glycoprotein of SARS-CoV-2 to appraise the impact of glycans in structure stabilization and antibody occlusion.

Structural and conformational analysis of SARS CoV 2 N-CTD revealing monomeric and dimeric active sites during the RNA-binding and stabilization: Insights towards potential inhibitors for N-CTD.

The advent of SARS-CoV-2 has become a universal health issue with no appropriate cure available to date. The coronavirus nucleocapsid (N) protein combines viral genomic RNA into a ribonucleoprotein and protects the viral genome from the host's nucleases. Structurally, the N protein comprises two independent domains: the N-terminal domain (NTD) for RNA-binding and C-terminal domain (CTD) involved in RNA-binding, protein dimerization, and nucleocapsid stabilization. The present study explains the structural aspects associated with the involvement of nucleocapsid C-terminal domain in the subunit assembly that helps the RNA binding and further stabilizing the virus assembly by protecting RNA from the hosts exonucleases degradation. The molecular dynamics (MD) simulations of the N-CTD and RNA complex suggests two active sites (site I: a monomer) and (site II: a dimer) with structural stability (RMSD: ~2 Å), Cα fluctuations (RMSF: ~3 Å) and strong protein-ligand interactions were estimated through the SiteMap module of Schrodinger. Virtual screening of 2456 FDA-approved drugs using structure-based docking identified top two leads distinctively against Site-I (monomer): Ceftaroline fosamil (MM-GBSA = -47.12 kcal/mol) and Cefoperazone (-45.84 kcal/mol); and against Site-II (dimer): Boceprevir, (an antiviral protease inhibitor, -106.78 kcal/mol) and Ceftaroline fosamil (-99.55 kcal/mol). The DCCM and PCA of drugs Ceftaroline fosamil (PC1+PC2 = 71.9%) and Boceprevir (PC1 +PC2 = 61.6%) show significant correlated residue motions which suggests highly induced conformational changes in the N-CTD dimer. Therefore, we propose N-CTD as a druggable target with two active binding sites (monomer and dimer) involved in specific RNA binding and stability. The RNA binding site with Ceftaroline fosamil binding can prevent viral assembly and can act as an antiviral for coronavirus.

Conformational flexibility and structural variability of SARS-CoV2 S protein.

Spike (S) glycoprotein of SARS-CoV2 exists chiefly in two conformations, open and closed. Most previous structural studies on S protein have been conducted at pH 8.0, but knowledge of the conformational propensities under both physiological and endosomal pH conditions is important to inform vaccine development. Our current study employed single-particle cryoelectron microscopy to visualize multiple states of open and closed conformations of S protein at physiological pH 7.4 and near-physiological pH 6.5 and pH 8.0. Propensities of open and closed conformations were found to differ with pH changes, whereby around 68% of S protein exists in open conformation at pH 7.4. Furthermore, we noticed a continuous movement in the N-terminal domain, receptor-binding domain (RBD), S2 domain, and stalk domain of S protein conformations at various pH values. Several key residues involving RBD-neutralizing epitopes are differentially exposed in each conformation. This study will assist in developing novel therapeutic measures against SARS-CoV2.

In Silico Structure-Based Repositioning of Approved Drugs for Spike Glycoprotein S2 Domain Fusion Peptide of SARS-CoV-2: Rationale from Molecular Dynamics and Binding Free Energy Calculations.

The membrane-anchored spike (S) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a pivotal role in directing the fusion of the virus particle mediated by the host cell receptor angiotensin-converting enzyme 2 (ACE-2). The fusion peptide region of the S protein S2 domain provides SARS-CoV-2 with the biological machinery needed for direct fusion to the host lipid membrane. In our present study, computer-aided drug design strategies were used for the identification of FDA-approved small molecules using the optimal structure of the S2 domain, which exhibits optimal interaction ratios, structural features, and energy variables, which were evaluated based on their performances in molecular docking, molecular dynamics simulations, molecular mechanics/generalized Born model and solvent accessibility binding free energy calculations of molecular dynamics trajectories, and statistical inferences. Among the 2,625 FDA-approved small molecules, chloramphenicol succinate, imipenem, and imidurea turned out to be the molecules that bound the best at the fusion peptide hydrophobic pocket. The principal interactions of the selected molecules suggest that the potential binding site at the fusion peptide region is centralized amid the Lys790, Thr791, Lys795, Asp808, and Gln872 residues.IMPORTANCE The present study provides the structural identification of the viable binding residues of the SARS-CoV-2 S2 fusion peptide region, which holds prime importance in the virus's host cell fusion and entry mechanism. The classical molecular mechanics simulations were set on values that mimic physiological standards for a good approximation of the dynamic behavior of selected drugs in biological systems. The drug molecules screened and analyzed here have relevant antiviral properties, which are reported here and which might hint toward their utilization in the coronavirus disease 2019 (COVID-19) pandemic owing to their attributes of binding to the fusion protein binding region shown in this study.